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ABSTRACT 

Context. Stellar activity induces apparent radial velocity (RV) variations in late-type main-sequence stars that may hamper the detec- 
tion of low-mass planets and the measurement of their mass. 

Aims. We use simultaneous measurements of the active planet host star HD 189733 with high-precision optical photometry by the 
MOST satellite and high-resolution spectra by SOPHIE. We apply on this unique dataset a spot model to predict the activity-induced 
RV variations and compare them with the observed ones. 

Methods. The model is based on the rotational modulation of the stellar flux. A maximum entropy regularization is applied to find a 
unique and stable solution for the distribution of the active regions versus stellar longitude. The RV variations are synthesized con- 
sidering the effects on the line profiles of the brightness perturbations due to dark spots and bright faculae and the reduction of the 
convective blueshifts in the active regions. 

Results. The synthesized RV time series shows a remarkably good agreement with the observed one although variations on timescales 
shorter than 4-5 days cannot be reproduced by our model. Persistent active longitudes are revealed by the spot modelling. They rotate 
with slightly different periods yielding a minimum relative amplitude of the differential rotation of ACl/ii = 0.23 ± 0. 10. Moreover, 
several active regions with an evolution timescale of 2 - 5 days and an area of 0. 1 - 0.3 percent of the stellar disc are detected. 
Conclusions. The method proves capable of reducing the power of the activity-induced RV variations by a factor from 2 to 10 at 
the rotation frequency and its harmonics up to the third. Thanks to the high-precision space-borne photometry delivered by CoRoT, 
Kepler, or later PLATO, it is possible to map the longitudinal distribution of active regions in late-type stars and apply the method 
presented in this paper to reduce remarkably the impact of stellar activity on their RV jitter allowing us to confirm the detection of 
low-mass planets or refine the measurement of their mass. 

Key words, stars: late-type - stars: activity - stars: rotation - techniques: radial velocity - planetary systems - stars: individual 
(HD 189733) 

1. Introduction ties of the cross-correlation function used to measure the RV 

itself, such as the mean line bisector shift an d/or its FWHM 

. ^ Several studies have recently addressed the perturbation of (^ g^ iMelo et al. 1 120071: iBoisse et al. 2009; P ont et al. ll20Hl) . 

^ the RV induced by stellar activity. This jitter may severely ^he main limitation of such techniques is the lack of a tight 

H hamper the detection of low-mass planets or the measure- coiTelation between the RV perturbation and those proxies for 

. >^ , ment of the mass of those found by the inethod of the ^^e average spectral line distortion, as demonstrated by, e.g., 

transits. ^ " ' ^/^/^ »» ^ i ^/^,/^_ _ i 



^ (12010D, | Meunieretal.| mm, and Boisse et al. (201 1 ). The other group consists of the techniques 

| Meunieretal.| (l2010b) have considered the Sun as a star find- ^^^^ g^pj^it the modulation of the activity-induced perturbation 
mg that the reduction of convective blueshifts in active regions, ^y the stellar rotation. Most of the power due to the activity- 
due to the quenching of turbulent motions by the photosphenc induced perturbation is indeed concentrated at the rotation fre- 
magnetic fields, has a remarkable effect mducing variations with qugncies and its first two or three harmonics. Therefore, it is pos- 
an amplitude up to 5 - 10 m s' . Pumusque et al. | mj> have ^j^le to filter it out leaving the signal due to the reflex orbital mo- 
extended those studies introducing a model for the impact of ^^^^ induced by a planetary companion, unless its orbital period 
dark spots on the RV variations of solar-hke stars based on the f^jj^ ^j^^g ^^ such frequen cies, and to fit simultaneously activity 
statistics of sunspot groups. ^^^j planetary signals (see iBoisse et al.ll201lL for details). This 
Several techniques have been proposed to reduce the impact approach has been recently applied, among others, to detect a 
of the RV perturbation induced by stellar activity. They can be second non-transiting planet around the active stai- CoRoT-7 and 
classified into two broad groups. Those in the first group exploit 
the correlation between the RV perturbation and other proper- 



A. F. Lanza et al.: Synthesizing the activity-induced RV variations of HD 189733 



to refine the measure ment of the mass of its transiting tellur ic 
planet CoRoT-7b (see lOueloz et al. ll2009l:lBoisse et al. II20TI . 

In addition to those methods, another approach to reduce the 
impact of stellar activity can be based on the models of the ac- 
tive region distribution derived from high-precision photometry 
such as that provided by the space-borne telescopes of MOST, 
CoRoT, or Kepler for their target stars. From their light curves, it 
is possible to derive the longitudes of the active regions and their 
area, tracing their evolution along successive rotations. This in- 
formation can be used to m odel the RV perturbation induced by 
stellar activity. [Lanza et aTl ( 12010 ) have applied this approach to 
the light curve of CoRoT-7 to predict the power spectrum of the 
RV perturbation induced by stellar activity and derive the signif- 
icance of the signals attributed to the telluric planets CoRoT-7b 
and CoRoT-7c. The main limitation of that study was the lack 
of simultaneity between the CoRoT light curve and the HARPS 
RV time series analysed to extract the reflex motions induced by 
the two planets. 

Here we further explore this approach by means of simulta- 
neous photometric and RV measurements of HD 189733 mak- 
ing use of the MO ST hght curve and the SOPHIE RV time 
series published bv iBoisse et al. I (l2009h . They span about one 
month in the summer of 2007, corresponding to two full rota- 
tions of the star with a period of P,ot ~ 12 days (iHenrv & Winn I 
l2008h . HD 189733 is accompanied by a Jupiter-mass transiting 
planet with an orbital period of 2.218 days. With an effective 
temperature of ~ 5000 K and an average chromospheric index 
Shk = 0.525, this Kl main-sequen ce star is remarkably active 
dBouchvetal. 1 120051: iMoutou et alj 12007: 'Bo isse et al.ll2009l) . 
Its radial velocity (hereafter RV) measurements are then per- 
turbed with variations up to ~ 20 - 30 m s ', as evident in the 
residuals of the RV best fit of the planetary orbit (iBoisse et al. I 
I2009h . 

Our spot modelling technique allows us to retrieve the lo- 
cation and the evolution of the photospheric active regions of 
HD 189733 finding also evidence of a remarkable ste llar differ- 
ential rotation, as recently found bv lFares et al. I (1201 Ot) by means 
of Doppler Imaging techniques. The information on the starspot 
longitude and evolution is used to synthesize the RV variations 
that can be compared with the simultaneous SOPHIE measure- 
ments to assess the advantages and limitations of the technique. 
We find that the modulation of the RV perturbation induced 
by stellar activity is reproduced remarkably well on timescales 
comparable with the stellar rotation period but it is not possi- 
ble to model the variations with a timescale of 1 - 2 days, i.e., 
much shorter than the stellar rotation period. This is an unavoid- 
able consequence of the fact that the rotational modulation of 
the optical flux is used to map starspots. On the other hand, the 
main advantage of this method is the remarkable reduction of the 
perturbations at frequencies comparable with the stellar rotation 
frequency and its low-order harmonics, that is an important im- 
provement over previously proposed approaches. 

The next section describes the observational dataset. Then, 
we define our model: 1) how active regions locations are derived 
from lightcurve; 2) how we compute the distortion of the line 
profiles due to the brightness perturbations owing to dark spots 
and bright faculae and the reduction of the convective blueshifts 
in the active regions; and 3) how the distributions of active re- 
gions derived from the spot model fitted to the light curve are 
used to synthesize RV variations. The fourth section gives the 
parameters chosen for our modelling of HD 189733. The results, 
given in Section 5, are discussed in the final section. 



2. Observations 

2.1. MOST photometry 

In 2007, HD 189733 was monitored by the MOST 
(Microvariability and Oscillations of STars) microsatellite 
for 30.4994 days covering ~ 2.5 rotation cycles of the star from 
HJD 2454298.55874 to 2454329.05818. MOST was launched 
in 2003 and was designed to observe nearby, relative bright 
stars through a single broadband (350-7 00 nm) filter using a 
15/1 7. 3-cm Rurnak-M atsutov telescope (IWalker et al. I I2003L 
Matthews et al. 1 12004|) . The photometry was collected in the 
so-called MOST's Direct Imaging mode, where a defocussed 
stellar image was recorded on a subraster of a single CCD. 
The data were reduced as described in Rowe et al. (2006). The 
stellar fluxes were extracted by aperture photometry and the 
raw instrumental light curve was decoiTelated against the sky 
background and also against the location of the centroid of the 
point-spread function on the CCD. The satellite had an 820 km 
altitude Sun-synchronous polar orbit of 101.4 minutes and the 
data were binned to this period. The transits of the known planet 
were removed from the MOST lightcurve. This led to a time 
series consisting of 414 points. As seen in the top panel of 
Fig. S HD 189733 exhibited in 2007 a ~ 1% variation in flux, 
which is much less than the ~ 3% variation observed in 2006 b y 
the same telescope dMiller-Ricci et al. II2008I: ICroll et al. Il2007h . 
Moreover, a decreasing long-term trend was seen that was not 
exhibited by the comparison stars, suggesting that the drift was 
real. 

We assume the maximum observed flux at 
HJD 2454300.31968 as the unit of measure for the flux 
and normalize the light curve to that value. The mean error of 
the binned points is 1.15 x 10 '^ in relative flux units. However, 
since MOST photometry is non-differential and can experience 
drifts on daily or weekly timescales, such a value should be 
considered as a lower limit to the photometric precision and 
doe s not ta ke into account systematic errors. For example, 
Croll et al. I (12007) pointed out that the stray light background 
can be modulated at a period of 1 day and its first harmonic (due 
to the Sun-synchronous orbit). The flux level corresponding 
to the star without spots is unknown, therefore we assume the 
maximum observed flux as the reference level for the spot 
modelling, i.e., we assume that it coiTesponds to the flux of the 
unspotted star in the MOST passband. We can also mention 
that the M-dwarf companion of HD 189733 dBakos et al. Il2006h 
does not contaminate appreciably the photometric signal. 



2.2. SOPHIE radial velocity measurements 

SOPHIE is a high-resolution spectrograph mounted on the 1.93- 
m telescope at the Observatoire de Haute -Provence, France 
(iPerruchot et al. llIOOSLlBouchv et al. ll2009bl) . The SOPHIE ob- 
servations of HD 189733 were conducted simultaneously with 
the MOST photometry. Usually two exposures per night were 
gathered during more than one month (from 12 July to 23 August 
2007) to sample both the short orbital period (~ 2.2 days) and the 
longer rotation period of the star (~ 12 days). The observations 
were conducted in the high-resolution mode (A/AA ~ 75000) 
and with an average signal-to-noise ratio (SNR) per pixel of 
about 80 at /I - 550 nm. The journal of the obs ervations and 
the instrument setup were reported in Sect. 2 of IBoisse et al. I 
(2009). 

For the present work, the original spectra were reduced again 
with a new pipeline that, among others, corrects for the effects 
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of the charge transfer inefficiency of the CCDs (iBouchv et al. I 
l2009ah and removes part of the spectrum that may contain 
telluric lines before applying the correlation mask. This gives 
slightly different weights to the red and blue parts of the spec- 
trum for the determination of the RV through the standard cross- 
correlation method, resulting in a RV offset of +21.6 m s~ ' with 
respect to the RV values reported in Boisse et al. (2009). The 
recently reported seeing effect (Boisse et al. 2010) was not well 
constrained and observed in these data, therefore it is not cor- 
rected. An average error of 3.9 m s ' for the RV measurements 
is then a good estimate of the accuracy of our values account- 
ing for the photon-noise error, the uncertainty in the wavelength 
calibration, and SOPHIE systematics. The effects of the system- 
atics is taken to be 3 m s ' for all the measurements, derived 
from the standard deviation of the RV of the most stable stars 
observed with SOPHIE, some of which have been observed with 
HARPS to have a stability better than 1 m s"'. These stars have 
more than 10 measurements and have been observed for more 
than six months in order to take into account the long-term sys- 
tematics. We note that new fiber scramblers have just been in- 
stalled on SOPHIE (June 2011), and preliminary tests showed 
a significant improvement in the RV accuracy. To compute the 
RV perturbations due to stellar activity, the orbital solution has 
been subtracted. Th e adopted orbital elements are those given by 
lBoisseetal.1(l2009h . 

3. Light curve and RV modelling 

3.1. Spot modelling of wide-band light curves 

The reconstruction of the surface brightness distribution from 
the rotational modulation of the stellar flux is an ill-posed prob- 
lem, because the variation of the flux vs. rotational phase con- 
tains information only on the distribution of the brightness inho- 
mogeneities vs. longitude. The integration over the stellar disc 
effectively cancels any latitudinal information, particularly when 
the inclination of the rotation axis along the line of sight is 
close t o 90°, as in the present case (see Sect. |4]and iLanza et aU 
l2009ah . Therefore, we need to include a priori information in 
the light curve inversion process to obtain a unique and stable 
map. This is done by computing a maximum entropy (hereinafter 
ME) map, which has been proven to best reproduce active re- 
gion distribution and area variations in the case of the Sun (cf. 
iLanza et ani2007h . For a di fiferent modell i ng app roach, based on 
discrete starspots, see, e.g..'Mosser et al.l (l2009h . 

In the present model, the stellar surface is subdivided into el- 
ements, i.e., in 200 squares of side 18°, with each element con- 
taining unperturbed photosphere, dark spots, and facular areas. 
The fraction of the k-th element covered by dark spots is indi- 
cated by its filling factor /;;, the fractional area of the faculae is 
Qfk, and the fractional area of the unperturbed photosphere is 
1 - {Q + l)fk. The contribution to the stellar flux coming from 
the k-th surface element at time tj, where j = 1 , ..., A^, is an index 
numbering the A^ points along the light curve, is given by: 



AFkj = /o(W;) {1 - (fi + l)fk + cjk+ 

Qfkli +Cf(l -tJ^kj)]}AkiUkjw(iJkj), 



(1) 



where Iq is the specific intensity in the continuum of the unper- 
turbed photosphere at the isophotal wavelength of the observa- 
tions, Cs and Cf are the spot and facular contrasts, respectively (cf. 
iLanza et ani2004l) . Ak is the area of the A:-th surface element. 



is its visibility, and 

fikj = cos i/r^y - sin / sin 6k cos [4 + ^(tj - fo)] + cos / cos Ok, (3) 

is the cosine of the angle i//kj between the normal to the surface 
element and the direction of the observer, with / being the incli- 
nation of the stellar rotation axis to the line of sight, 6k the colat- 
itude and (k the longitude of the k-th surface element, Q. the an- 
gular velocity of rotation of the star (O = In /Prod, and to the ini- 
tial time. The specific intensity in the continuum va ries accord- 
ing to a quadratic limb-darkening law, as adopted bv .Lanza et al.l 
(l2003h for the case of the Sun, viz. /o oc Op -h bpfi H- c^fi^. The stel- 
lar flux computed at time tj is then: F{tj) = Y^k ^^kj- To warrant 
a relative precision of the order of 1 0"^^ in the computation of the 
flux F, each surface element is further subdivided into elements 
of 1 ° X 1 ° and their contributions, calculated according to Eq. ([T]i, 
are summed up at each given time to compute the contribution 
of the 18° X 18° surface element to which they belong. 

We fit the light curve by changing the value of the spot fill- 
ing factor / over the surface of the star, while Q is held constant. 
Even fixing the rotati on period, the incli nation, and the spot and 
facular contrasts (see iLanza et aDl2007l for details), the model 
has 200 free parameters and suffers from non-uniqueness and 
instability. To find a unique a nd stable spo t map, we apply ME 
regularization, as described in lLanzaet alJ (l2007h . by minimiz- 
ing a functional Z, which is a linear combination of the x^ and 
the entropy functional S ; i.e., 

z^x\f)-nS(f), (4) 

where / is the vector of the filling factors of the surface ele- 
ments, 7/ > a Lagrangian multiplier determining the trade-off 
between light curve fitting and regularization, and the expres- 
sion of 5 is: 



z 



Wk 



m 1 - m 



(5) 



w(Mkj) 



1 ifiukj>0 

if;U*;<0 



(2) 



where Wk is the relative area of the k-th surface element (total 
surface area of the star - 1) and m the default spot covering 
factor that fixes the limiting values for fk as: m < fk < ( 1 - m). In 
our modelling we adopt m - 10"^ (cf. iLanza et al] |l998^. The 
entropy functional S is constructed in such a way that it attains 
its maximum value when the star is immaculate. Therefore, by 
increasing the Lagrangian multiplier t], we increase the weight of 
S in the model and the area of the spots is progressively reduced. 
This gives rise to systematically negative residuals between the 
observations and the best-fit model when 77 > 0. 

The optimal value of rj depends on the information content 
of the light curve, which in turn depends on the ratio of the am- 
plitude of its rotational modulation to the average standard de- 
viation of its points. In the case of HD 189733, the amplitude 
of the rotational modulation is ~ 0.015, while the nominal stan- 
dard deviation of the points is ~ 1.15 x 10 '^ in relative flux units 
(see Sect. 12. lb . giving a signal-to-noise ratio of ~ 130. However, 
since short-term instrumental fluctuations cannot be excluded in 
the case of MOST photometry (cf. Sect. 12.1b . we prefer to con- 
sider the standard deviation of the residuals obtained with the 
unregularized best fit (i.e., with 77 = 0; see Sect. 15.1b . In such a 
way, we find that the relative accuracy of MOST photometry is 
overestimated by at least a factor of ~ 3, i.e. the mean standard 
deviation of the photometric points is ~ 4 x 10 '*, leading to a 
signal-to-noise ratio of ~ 40 - 50. 

To fix the optimal value of the Lagrangian multiplier 77, we 
compare the modulus of the mean of the residuals of the regular- 
ized best fit |//regl with the standard error of the residuals them- 
selves, i.e., eo = ctq/ Va^, where ctq is the standard deviation of 
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the residuals of the unregularized best fit and A^ is the number 
of points in each fitted subset of the light curve having a dura- 
tion Aff (see below). We iterate until \i-ireg\ - jSfo, where y6 is a 
numerical parameter that will be fixed a posteriori according to 
the requisites of an acceptable quality for the fit and a regular 
evolution of the spot pa ttern from one spot map to the next (cf . 
the case of CoRoT-4 in iLanza et al. I l2009b'). 

In the case of the Sun, by assuming a fixed distribution of the 
filling factor, it is possible to obtain a good fit of the irradiance 
changes only for a limited time interval Aff, not exceeding 14 
days, which is the lifetime of the la rgest sunspo t groups domi- 
nating the irradiance variation ( Lanza et al.ll2003h . In the case of 
other active stars, the value of Aff must be determined from the 
observations themselves, looking for the maximum data exten- 
sion that allows us a good fit with the applied model (see Sect.]?]). 

The optimal values of the spot and facular contrasts and of 
the facular-to-spotted area ratio Q in stellar active regions are 
unknown a priori. In our model the facular contrast Cf and the 
parameter Q enter as the product CfQ, so we can fix C{ and vary 
Q, estimating its best value by minimizing the x^ of the model, 
as shown in Sect. |4] Since the number of free para meters of the 
ME m odel is large, we make use of the model of JLanzaet alJ 
(l2003h to fix the value of Q. It fits the light curve by assuming 
only three active regions to model the rotational modulation of 
the flux plus a uniformly distributed background to account for 
the variations of the mean light level. This procedure is the same 
as adopted for CoRoT-2 and CoRoT-4 to fix the value of Q (cf. 
ILanza et al.ll2009aHLanza et al. Il2009bh . 

We shall assume an inclination of the rotation axis of 
HD 189733 of / ^ 85° (see Sect. Hi. Since the information on 
spot latitudes that can be extracted from the rotational modula- 
tion of the flux for such a high inclination is negligible, the ME 
regularization virtually puts all the spots at the sub-observer lat- 
itude (i.e., 90° - / ~ 5°) to minimize their area and maximize the 
entropy. Therefore, we are limited to mapping only the distribu- 
tion of the active regions vs. lon gitude, which can be done with 
a resolution of at least ~ 50° (cf. Lanza et"an i2"007. 2009a). Our 
ignorance of the true facular contribution to the light modulation 
may lead to systematic errors in the a ctive region longitu des de- 
rived by our model, as discussed by iLanza et al.l (l2007h in the 
case of the Sun. 



3.2. Modelling the line distortion induced by activity 

Surface magnetic fields in late-type stars produce bright- 
ness and convection inhomogeneities that distort their spec- 
tral line profiles leading to apparent R V variations (cf., e.g., 
ISaar & Donahuelll997l;[Saar et al.|[T998l) . To compute the appar- 
ent RV variations induced by stellar active regions, we adopt a 
simple model for each line profile considering the residual pro- 
file R{A) at a wavelength A along the line, i.e.: R{A) = [1 - 
/(/l)//c], where I{A) is the specific intensity at wavelength A and 
/c the intensity in the continuum adjacent to the line. The local 
residual profile is assumed to be a Gaussian with thermal and 
microturbulent width A/Id (cf., e.g., iGravii il988i, Ch. 1), i.e.: 



R(A) oc exp 



A-Aq 
AAn 



(6) 



where Aq is the local central wavelength. The different effective 
temperatures in dark spots or bright faculae modify the depth of 
the local line profile. To take into account this effect, we mul- 
tiply the unperturbed profile given by Eq. (|6]l by (1 - Cifu), 




0.985 



O.O 0.2 0.4- 0.6 0.8 1 .O 
Rotation phase 

Fig. 1. Upper panel: The RV variation produced by a single 
equatorial spot on HD 189733 vs. the rotation phase for different 
values of the parameter Cd according to the different linestyles: 
Cd - 0.0 (no temperature effect on the line profile): solid line; 
Cd - 0.2: dotted line; and Cd = 0.4: dashed line. Lower panel: 
the corresponding modulation of the optical flux vs. the rotation 
phase. Note that the three plots overlap because the modulation 
is independent of Cd (see the text). 



where Cd is a coefficient that depends on the specific line con- 
sidered and the physical properties of the stellar atmosphere, and 
fk is the spot filling factor in the specified surface element. This 
factor accounts for the temperature reduction in dark spots, but 
does not consider the increased effective temperature in facu- 
lae because the facular contrast is assumed to be solar-like (cf. 
Sect. H, i.e., corresponding to a modest temperature increase 
of » 100 K in the layers where the optical spectrum is pro- 
duced (cf., e.g., lUnruh et al. II 199 9V Since the RV measurement 
is based on a cross-correlation function, Cd must be averaged 
over a great number of spectral lines having in general quite dif- 
ferent temperature responses. For this reason, we consider Cd as 
a free parameter in our model and adjust it to fit the observed 
RV variations induced by stellar activity (see below for the ef- 
fect of the variation of Cd on the RV perturbation induced by a 
starspot). In addition to the variation of the local line profile, the 
local specific intensity along a spectral line I{A) depends on the 
variation of the continuum intensity /c owing to limb-darkening 
and the effects of dark spots and bright faculae, as specified by 
Eq. ©. 

To include the effects of surface magnetic fields on convec- 
tive motions, we consider the decrease of macroturbulence ve- 
locity and the reduction of convective blueshifts in active re- 
gions. Specifically, we assum e a lo cal radial-tangential macro- 
turbulence, as introduced by iGravl ( Il988.) . with a distribution 
function for the radial velocity v of the form: 



(7) 




®{v,iD = - 



^Jn^KTU 



where ^rj is the macroturbulence and fi is given by Eq. (O. 
The profile emerging from each surface element is obtained by 
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convolving the local profile in Eq. (|6]l with the Doppler shift 
distribution as generated by the macroturbulence function in 
Eq. Q. The parameter ^rt is assumed to be reduced in spot- 
ted and facular areas, according to their total filling factor, i.e., 
fc - fo[l - (2 + 1)/]. where (o is the unperturbed macroturbu- 
lence. This is a convenient parameterization of the quenching ef- 
fect of surface magnetic fields on convective turbulence, at least 
in the framework of our simplified model, and is supported by 
the observations of the reduction of t urbulent velocity fields in 
the plage regions of the Sun (cf., e.g., iTitle et aDll992h . 

Convective blueshifts arise because in stellar photospheres 
most of the flux comes from the extended updrafts at the cen- 
tre of convective granules. At the centre of the stellar disc, the 
vertical component of the convective velocity produces a max- 
imum blueshift, while the effect vanishes at the limb where the 
projected velocity is zero. The cores of the deepest spectral lines 
form in the upper layers of the photosphere where the vertical 
convective velocity is low, while the cores of shallow lines form 
in deeper layers with higher vertical velocity. Therefore, the 
cores of shallow lines are blueshifted with respect to the cores 
of the deepest lines. Gray (2009) shows this effect by plotting 
the endpoints of line bisectors of shallow and deep lines on the 
same velocity scale. He shows that the amplitude of the relative 
blueshifts scales with the spectral type and the luminosity class 
of the star For the G8V star r Ceti, which has a spectral type 
not too much different from the KIV of HD 189733, the con- 
vective blueshifts are approximately similar to those of the Sun, 
so we adopt solar values in our simulations. In active regions, 
vertical convective motions are quenched, so we observe an ap- 
parent redshift of the spectral lines in spotted a nd facular areas 
in comp arison to the unperturbed photosphere. iMeunier et al. I 
(1201 Obi) quantify this effect in the case of the Sun, and we adopt 
their results, considering an apparent redshift AVf - 200 m s"' 
in faculae and AVs - 300 m s"' in cool spots. 

In principle, the integrated effect of convective redshifts can 
be measured in a star by comparing RV measurements obtained 
with two different line masks, one in cluding the shallow lines 
and the other the deep lines (cf. Meuni er et al. 1 1201 Ob). For 
HD 18 9733, which lacks such measurements, we apply the re- 
sults of lGravl ( l2009h and adopt solar-like values as the best ap- 
proximation. 

Considering solar convection as a template, intense down- 
drafts are localized in the dark lanes at the boundaries of the 
upwelling granules, but they contribute a significantly smaller 
flux because of their lower brightness and smaller area. While 
a consideration of those downdrafts is needed to simulate the 
shapes of line bisectors, it is beyond the scope of our simplified 
approach that assumes that the whole profile of our template line 
forms at the same depth within the photosphere. Therefore, we 
restrict our model to the mean apparent RV variations by ne- 
glecting the associated variations of the bisector shape and do 
not include the effect of convective downdrafts, as well as those 
of other surface flows typical of solar active regions, such as the 
Evershed effect in sunspots (cf., e.g., iMeunier et al. 11201 Ob*) . 

The local central wavelength Aq of the ^-th surface element 
at time tj is given by Ao^j = /Ir(1 + Vkjic), where v^j is its radial 
velocity. 



vective blueshifts in spots and faculae, which is parameterized 

as 



"^kj - -{v sin sin Ok sin[4 + ^(tj - to)] + AVcs, 



(8) 



AVcs = (Av, + Q^v^fk^ikj, 



(9) 



where fk is the spot filling factor of the element, Q the facular- 
to-spotted area ratio, and Hkj is given by Eq. ([3]). 

We integrate the line specific intensity at a given wavelength 
over the disc of the star using a subdivision into 1 ° x 1 ° surface el- 
ements to obtain the flux along the line profile F(A, tj) at a given 
time tj. To find the apparent stellar RV, we can fit a Gaussian to 
the line profile F{A, tj), or we can determine the centroid of the 
profile as 



^ce(f;) = 



^ A'R{A,tj)dA 

f'R{A,tj)dA ' 



(10) 



with c the speed of light, v sin / the stellar projected rotational 
velocity, Ar the rest wavelength of the spectral line, and AVcs the 
apparent convective redshift arising from the reduction of con- 



where 'R= [I -F(A, tj)IFJ^tj)] is the residual profile of the spec- 
tral line computed from the ratio of the flux in the line F{A, tj) to 
the flux /^c(0) ill the adjacent continuum at any given time tj. For 
HD 189733, the two methods lead to almost equal values of the 
RV perturbation with differences never exceeding ~ 3 percent. 

A single line profile computed with the above model can 
be regarded as a cross-correlation function (hereafter CCF) ob- 
tained by cross-correlating the whole stellar spectrum with a line 
mask consisting of Dirac delta functions giving t he rest wave- 
length and depth of each individual line (e.g., lOueloz et al. I 
1200 Ih . Therefore, to derive the RV from a single synthetic Une 
profile is equivalent to measuring the RV from a CCF, either 
by fitting it with a Gaussian or by computing its centroid wave- 
length. A better approach would be to simulate the whole stel- 
lar optical spectrum or an extended section of it to account for 
the wavelength dependence of the spot and facular contrasts, as 
well as of the convective inhomogeneities (cf., e.g., 'Deso rt et al.l 
2007; Lagrange et al. 2010; Meunier et al. 2010b). Again, in 
view of our limited knowledge of the distribution of active re- 
gions on the stellar surface, we believe that our simplified ap- 
proach is adequate for estimating the RV perturbations induced 
by magnetic fields. 

An illustration of the effects of dark spots, bright faculae, and 
macroturbulen ce perturbations on the stellar RV measurements 
is presented in iLanza et aP (12010 ) together with a comparison 
of the results of the present model with those of previous ap- 
proaches. Here we only show in Fig. [1] (upper panel) the effect 
of the variation of the paramet er Ca on the RV cur ve of a single 
spot that was not considered bv lLanza et al.l ( l2010l) . The adopted 
spot is located on the equator of the star and has a square shape 
with a side of 18°, a covering factor / - 0.99, and a contrast 
Cs = 0.482 without any facular contribution. The stellar and line 
profile parameters are assumed to be the same as adopted for 
the modefling of the RV variations of HD 189733 (see Sect.©. 
The increase of the amplitude of the RV variation with Cd is a 
consequence of the higher bump associated with the spot on the 
line profile when the depth of the local profile is reduced by the 
factor (1 - Cd/). On the other hand, the modulation of the opti- 
cal flux associated with the transit of the spot across the stellar 
disc is independent of the value of Cd because it depends on the 
flux in the continuum outside the spectral line (cf. Fig.[Tl bottom 
panel). 

The amplitude of the RV variation produced by an active 
region depends on several parameters that are poorly known, 
i.e., the latitude of the active region, the spot and facular con- 
trasts, and the macroturbulence parameter which is difficult to 
disentangle from rotati onal broadening in a slowly rotating star 
such as HD 189733 dLeger et al. II2009I) . Moreover, the spot and 
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facular contrasts depend on the wavelength (iLanza et al.ll2004l) . 
leading to a difference of « 10 percent in the RV variations as 
derived from diffe rent orders of an echelle spectrum (cf., e.g., 
iDesort et al.ll2007h . The simultaneous presence of several ac- 
tive regions gives rise to a complex line profile distortion in the 
case of a slowly rotating star because the perturbation of each 
active region is not separated from the other by the rotational 
broadening. This implies an additional 10-15 percent uncer- 
tainty in the determination of the central wavelength of the pro- 
file by the Gaussian fit or the centroid method, even for an in- 
finite signal-to-noise ratio and perfectly regular sampling, as in 
the case of simulated line profiles. In consideration of all these 
limitations, the absolute values of the RV variations computed 
with our model are uncertain by 10-20 percent in the case of 
complex distributions of active regions, such as those derived by 
our spot modelling technique as applied to MOST and CoRoT 
light curves. 

3.3. RV variations from spot modelling 

We can use the distribution of active regions as derived from our 
spot modelling to synthesize the corresponding RV variations 
according to the approach outlined in Sect. 13.21 Since our spot 
model assumes that active regions are stable for the time interval 
of each fitted time series Aft , the distribution of surface inhomo- 
geneities can be used to synthesize the RV variations having a 
timescale of Aff or longer. To improve the time resolution, we 
shall consider models obtained from time intervals whose ini- 
tial time is shifted by Aff /2 from the previous interval, and av- 
erage the RV variations computed at the same time to account 
for the spot evolution in a simple way. Our ME modelling mini- 
mizes the area of the starspots and puts them close to the equator. 
Since the spot latitudes cannot be determined from photometry 
(cf. Sect. 13.1b . we assume that all the spots are located close to 
the equator, which maximizes the computed RV variations (cf 
Eq.[8]i. Note that the frequency content of the simulated RV time 
series is not significantly modified by this assumption because 
the duration of the transit of a given spot across the stellar disc is 
pratically independent of its latitude given that the inclination of 
the stellar rotation axis is close to 90° and our reference frame is 
rigidly rotating with the mean stellar rotation period. Moreover, 
the occurrence times of the maximum and minimum of the RV 
variations associated with a given spot are mainly determined by 
its longitude and its projected area that our modelling approach 
consistently retrieves from the modulation of the optical flux. 
Therefore, we do not expect any significant modification of the 
power spectrum of the synthetic RV variations associated with 
the assumption of equatorial spots. 

The active regions with lifetimes shorter than Aff /2 produce 
a photometric modulation that appears in the residuals of the 
best fit to the lightcurve. As we shall see in Sect. 15. II most of the 
short-term variability occurs on timescales of 1 -2 days, i.e., sig- 
nificantly shorter than the rotation period of HD 189733, so we 
can neglect, as a first approximation, the variation in the disc po- 
sition of active regions due to stellar rotation and estimate their 
area as if they were located at the centre of the disc. Of course, 
the actual RV perturbation depends on the position of the active 
regions on the disc, but this information is lacking. Therefore, 
the approach we describe below cannot be used to synthesize 
the actual RV variation, but only to estimate its maximum am- 
plitude. 

First, we express the photometric residuals as the relative de- 
viation AF/Fo between the observed flux and its best fit mea- 
sured in units of the reference unspotted flux Fq. Then we sub- 



tract its mean value fi^^s = (AF/Fq) from AF/Fq because the 
mean value corresponds to a uniformly distributed population of 
active regions that does not produce any RV variation. At each 
given time f, we adopt \AF/Fo - jUreskjT' as a measure of the 
filling factor of the active regions producing the short-term RV 
variations. In doing so, we neglect limb-darkening effects and 
assume that those active regions consist of dark spots with a 
contrast Cj. Finally, we compute the amplitude of the RV per- 
turbation due to such brightn ess inhomogeneities by means of 
Eq. (1) of lDesortetaP (12007 ) obtaining 



AV, 



Ri 



800 



AF 



c:\vsmiy\ 



(11) 



where the RV perturbation is measured in m s ' and the vsin/ 
of the star in km s"' . Eq. ( fTTT ) gives an upper limit for the am- 
plitude of the RV perturbation because it assumes that an active 
region spans the diameter of the stellar disc. On the other hand, 
the effects of convective redshifts in short-lived active regions is 
estimated as 



AVrc = (AV, + QAVf) 



AF 



-yUn 



(12) 



By adding both short-term perturbations, we compute the ampli- 
tude of the total variation induced by the active regions evolving 
on time scales shorter than the time resolution of our model, and 
add it to the standard deviation of the RV variations computed 
with the slowly evolving active regions to compute the total un- 
certainty of the synthesized RV values. 

4. Model parameters 



The basic stellar parameters are taken from Bouchv et al. ("20051) 
and Triaudetal. (2009) and are listed in Table [1] The limb- 
darkening parameters have been derived from iKuruczl (120001) 
model atmospheres for Teff = 5050 K, \ogg = 4.53 (cm s"^) and 
solar abundances, by adopting the MOST transmission profile in 
IWalker et al. I (12003 ). We compared our limb-darkening profile 
with that derived by Miller-Ricci et al. (2008), who adopted a 
different limb-darkening law and used MOST data acquired in 
2006, and found a maximum relative difference of 3 percent. 

The rotation period adopted for our spot modelling has been 
de rived fr om the optical light curve (Piot = 11 .953 ± 0.009 days) 
bv lHenrv & Winn (2008) . The polar flattening of the star due to 
the centrifugal potential is computed in the Roche approxima- 
tion with a rotation period of 1 1.95 days. The relative difference 
between the equatorial and the polar radii is e = 2.58 x 10"^ 
which induces a completely negligible relative flux variation of 
» 10"^ for a spot coverage of ~ 2 percent, as a consequence of 
the gravity darkening of the equatorial regions of the star 

The inclination of the stellar rotation axis is constrained 
by the observation s of the Rossiter-McLaughlin effect by 
iTriaud et al. I (12009 ^) who found that the angle y6RM between the 
projected directions on the plane of the sky of the orbital angular 
momentum and the stellar spin is j3rm - V-^5 + 0?32. Although 
the possibility that the stellar spin be not aligned with the orbital 
angular momentum cannot be excluded with certainty because 
only the sky-projected angle between the two vectors has been 
measured, we regard this occurrence as highly improbable and 
assume that the rotation axis is inclined with respect to the line 
of sight of the same angle as the orbit normal. 

The maximum time interval that our model can accurately 
fit with a fixed distribution of active regions Aff has been de- 
termined by dividing the total interval, T - 30.4994 days, into 



A. F. Lanza et al.: Synthesizing the activity-induced RV variations of HD 189733 



Table 1. Parameters adopted for the light curve and radial veloc- 
ity modelhng of HD 189733. 



Parameter 




Value 


Ref." 


Star Mass (Mq) 




0.823 


T09 


Star Radius (Ra) 




0.766 


T09 


Teff (K) 




5050 


B05 


log^(cms"^) 




4.53 


BOS 


Qp 




0.213 


LU 


bp 




0.933 


Lll 


Cp 




-0.152 


Lll 


/>,,, (days) 




11.953 


HW08 


e 


8.54 


X 10-* 


Lll 


Inclination (deg) 




85.3 


B05, T09 


Cs 




0.482 


Lll 


Ct 




0.115 


L04 


Q 




0.0 


Lll 


Atf (days) 




7.625 


Lll 


vsin; (km s"') 




3.316 


T09 


AiD(kms-i) 




2.33 


Lll 


fo (km s"') 




1.2 


Lll 




jferences: B05:'Bouchv et al. ' 


■200^: HW08:lHenrv& Winn 1 


(I2008,);L04: .Lanza etal. 


(200' 


l);Lll:i: 
U2OO90. 


)resent study; T09: 


iTriauc 


etal 





A^f equal segments, i.e., Aff - T/Nf, and looking for the mini- 
mum value of Nf that allows us a good fit of the light curve, as 
measured by the x^ statistics. We found that for Af < 4 the qual- 
ity of the best fit degrades significantly with respect to higher 
values, owing to a substantial evolution of the pattern of sur- 
face brightness inhomogeneities. We adopt At = 4, and therefore 
Aff - 7.625 days is the maximum time interval to be fitted with 
a fixed distribution of surface active regions in order to estimate 
the best value of the parameter Q (see below). 

To compute the spot contrast, we adopt the same mean tem- 
perature difference as derived for sunspot groups f rom their 
bolometric contrast, i.e. 560 K (iChapman et al.lll994l) . In other 
words, we assume a spot effective temperature of 4490 K, 
yielding a contrast Cs = 0.482 in the MOST passband (cf. 
[Lanza et al. 2007). A different spot contrast changes the absolute 
spot coverages, but does not significantly affect thei r longitudes 
and the ir time evolution, as discussed in detail bv iLanza et al.l 
(l2009ah . Therefore, adopting a spot temperature deficit of, e.g., 
~ 1000 K, as discuss ed in the analysis of HD 189733 HST data 
bv lPont et al. I (l2007h . will decrease the filling factor in propor- 
tion to the lower contrast, but does not change our results in any 
significant way. 

In the models computed with a facular contribution, the 

facular contr ast is assumed to be solar-like with Cf - 0.115 

dLanza et al.ll2004.) . The best value of the area ratio Q between 

the faculae and the spots in each active re gion has been esti- 

mated by means of the three-spot model bv iLanza et al.l (l2003l 

cf. Sect. 13.1b . In Fig. |2] we plot the ratio x^/x^ of the total x^ 

of the composite best fit of the entire time series to its minimum 

value y^ . , versus Q, and indicate the 95 percent confidence level 

I — II 1 

as derived from the F-statistics (e.g., iLampton et al.lll976l) . The 

choice of Aff = 7.625 days allows us to fit the rotational mod- 
ulation of the active regions for the longest time interval during 
which they remain stable, modelling both the flux increase due 
to the facular component when an active region is close to the 




Fig. 2. The ratio of the x^ of the composite best fit of the entire 
time series of HD 1 89733 to its minimum value vs. the parameter 
Q, i.e., the ratio of the area of the faculae to that of the cool 
spots in active regions. The horizontal dashed line indicates the 
95 percent confidence level for;tf^/;t^^ , determining the interval 
of acceptable Q values. 



limb as well as the flux decrease due to the dark spots when the 
same region transits across the central meridian of the disc. In 
such a way, a measure of the relative facular and spot contribu- 
tions can be obtained leading to an estimate of Q. Fig. |2] shows 
that the best value of Q is Q = 0, with an acceptable range ex- 
tending from ~ to ~ 5. Therefore, we adopt Q - for all our 
modelling in Sect. |5]unless otherwise indicated. We shall further 
comment on the value of Q in Sect. |6] 

To compute the RV variations induced by surface inhomo- 
geneities, we assume a line rest wavelength of 600 nm and a lo- 
cal thermal plus microturbulence broadening AAu - 2.3 km s"' . 
The vsin/ = 3.316 + 0.07 km s"' is deri ved from the obser - 
vations of the Rossiter-McLaughlin effect dTriaud et al. 1120091) . 
Although such a measurement can be affected by a system- 
atic error, we prefer to use this value, that is l arger than the 
2.97 + 0.22 km s' adopted by Winn et al] (120061) . to maximize 
the synthesized RV variations. Adopting vsin; = 2.97 km s' 
would imply an amplitude of the RV variations lower by a; 10 
percent for the same value of the parameter d. On the other 
hand, we could still reproduce the observed variations by chang- 
ing the average value of Cd (see Sect. I5.3l l. so that the choice of 
the V sin / value is not very critical for our modelling. The radial- 
tangential inacroturbulence velocity is assumed ^o = 12 km s' 
after iGravl (Il988h because it is difficult to disentangle from the 
rotational broadening owing to the slow rotation of the star. 



5. Results 

5. 1 . Light curve models 

We applied the model of Sect. l3.1l to the out-of-transit light curve 
of HD 189733, considering time intervals Aff - 7.625 days. The 
initial epoch of each fitted interval is shifted by Aff/2 = 3.81 
days, to sample the evolution of the starspots in a more con- 
tinuous way. The best fit without regularization (rj - 0) has a 
mean yUies = 2.15 x 10"^ and a standard deviation of the resid- 
uals (To = 3.99 X lO"'* in relative flux units. The Lagrangian 
multiplier 77 is iteratively adjusted until the mean of the resid- 
uals //res = -7.82 X lO-'^ - 2cro/ Va^, where N = 103 is the 
mean number of points in each fitted light curve interval Aff; the 
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standard deviation of the residuals of the regularized best fit is 
cr - 4.86 X 10""*. In other words, we fixed/? ^ 2 in the regular- 
ization procedure which gives the best compromise between the 
accuracy of the fit and the smoothness of the spot distribution 
and evolution. 

The composite best fit to the entire light curve is shown in the 
upper panel of Fig. [3] while the residuals are plotted in the lower 
panel. The computed flux values at the same time from succes- 
sive spot models are averaged. In spite of that, at the matching 
points between successive best fits we observe some disconti- 
nuities in the first derivative of the flux variation vs. time. This 
may in principle affect the rate of change of the RV perturba- 
tion induced by the spots vs. the rotation phase. However, we do 
not worry about correcting for this effect because the discontinu- 
ities always occur at epochs when no RV observation has been 
obtained. Therefore, the synthesizes RV values to be compared 
with the SOPHIE RV time series are not appreciably affected by 
those discontinuities. 

The residuals show oscillations with a typical timescale of 
~ 1 - 2 days that can be related to the rise and decay of active 
regions that cover » 0.2 - 0.3 percent of the stellar disc, i.e., 
comparable with the largest sunspot groups. The projected area 
of those active regions is estimated from the amplitude of the 
flux residuals and the adopted spot contrast, while their lifetimes 
are estimated by the duration of the residual fluctuations. These 
active regions cannot be modelled by our approach because they 
do not produce a significant rotational flux modulation during 
the ~ 12 days of the stellar rotation period as they move across 
the disc by only « 60° in longitude. Their estimated size and 
lifetime are not significantly degenerate with each other because 
our model has 200 degrees of freedom to fit any minute flux vari- 
ation with a timescale of several days by adjusting the covering 
factors /t of its 18° x 18° surface elements. In other words, by 
subtracting our model we filter out all the variations with time 
scales longer than 4-5 days leaving only the effects of the ac- 
tive regions evolving on shorter timescales. 

By decreasing the degree of regularization, i.e., the value of 
;;, we can marginally improve the best fit, but at the cost of in- 
troducing several small active regions that wax and wane from 
one Aff time interval to the next and are badly constrained by the 
rotational modulation. Nevertheless, the oscillations of the resid- 
uals do not disappear completely even for ?/ = 0, indicating that 
HD 189733 has a population of short-lived active regions with 
typical lifetimes of 1 - 2 days. Spots can be mapped along the 
strip occulted by the planet, as shown by Pont et al. I (l2007l) . but 
this is beyond the scope of the present investigation. 



5.2. Longitude distribution of active regions and stellar 
differential rotation 

The distributions of the spotted area vs. longitude are plotted in 
Fig-Sfor the seven mean epochs of our individual subsets adopt- 
ing a rotation period of 11.953 days. The longitude zero corre- 
sponds to the point intercepted on the stellar photosphere by the 
fine of sight to the centre of the star at HJD 2454298.55874, 
i.e., the sub-observer point at the initial epoch. The longitude in- 
creases in the same direction as the stellar rotation. This is con- 
sistent with the reference frames adopted in our previous studies 
JLanzaetal .11200931; iLanza et aril2009bh . but does not aflow a 
direct comparison of the mapped active regions with the dips in 
the light curve. 

Our maps show that the distribution of the spotted area ver- 
sus longitude evolves rapidly, which makes difficult to trace the 
evolution and migration of the active regions in an unambiguous 



way. Nevertheless, five main active regions can be identified in 
Fig. |4] and their migration is traced with dashed lines labelled 
with letters from A to E. Their rotation rates with respect to the 
adopted reference frame have been estimated by a linear best 
fit, i.e., by assuming a constant migration rate, and range from 
-5.9 + 1.5 deg/day for that labelled with D to 1.1 ± 1.5 deg/day 
for that labelled with A. Note that active longitudes D and E are 
not detectable in some of the ME distributions plotted in Fig.|4] 
but we decided to trace their evolution in a continuous way at- 
tributing the change of their visibility to a combination of their 
fast intrinsic evolution and the differential rotation. 

The relative amplitude of the surface differential rotation, 
as estimated from the difference between the maximum and the 
minimum migration rates reported above, is AQ/Q = 0.23+0.10. 
Unfortunately, no information is available on the spot latitudes, 
thus our AO/Q is only a lower limit. Since HD 189733 is more 
active than the Sun, in principle its active regions may cover a 
latitude range larger than in the Sun where sunspot gro ups ar e 
confined to ±40° from the equator (see, e.g., Strassmeier 20091) . 
iFares et alTI (|2010) using Zeeman Doppler Imaging and adopt- 
ing a solar-like differential rotation profile, estimated equatorial 
and polar rotation periods of 1 1 .94 ± 0. 1 6 and 1 6 .53 + 2 .43 days, 
giving a relative amplitude AQ/f2 = 0.38 + 0.18. Adopting their 
differential rotation profile with the above uncertainties, we es- 
timate a latitude ranging from 40° to 80° for the highest lati- 
tude spots, i.e., those in the active longitudes B and D. On the 
other hand, spots in the active longitudes A and C are close to 
the equator and those in E are at intermediate latitudes. Since 
the RV variations synthesized with all the spots located close to 
the equator reproduce the observations well (cf. Sect. 15.31 ). the 
lower value for the latitudes of the spots in B, D, and E, i.e., 
~ 30° - 40°, is favoured. Note that assuming, in the next sub- 
section, those spots close to the equator leads to an overestimate 
of their effects in the induced RV variations of ~ 20 - 25 percent 
which is within our typical errorbars of ±5 m s"' . 

The spot longitudes in our model show in general a de- 
pendence o n the facular-to-spotted area ratio Q (cf., e.g., 
Lanza et al] |2007h . therefore we computed regularized models 
with the maximum allowed facular contribution corresponding 
to 2 = 5. The relative amplitude of the differential rotation is 
not significantly modified with respect to the case wi t h Q = 0, 
i.e., it is still compatible with the result of lFares et al. I (l2010l) . 

The active longitudes traced in Fig. |4] show significant 
changes on a timescale as short as 3 - 4 days although their 
overall evolution may exceed the duration of the observations, 
i.e., ~ 30 days. The variation of the total spotted area vs. time 
is plotted in Fig. |5] and shows a monotonous increase along the 
MOST observing run corresponding to the overall decrease of 
the mean flux during the same interval. Note that the area per 
longitude bin and the total area depend on the spot and facular 
contrasts adopted for the modelling. Specifically, darker spots 
lead to a smaller total area while the effect of the facular contrast 
is more subtle and influences somewhat the longitudinal distri- 
bution of the active regions (see, Lanza et al. 2007, 2009ai, for 
more details). 

Finally, we note that the time resolution of our spot models 
is not adequate to look for a possible modulation of the stellar 
activity with the orbital period of it s close-in massive planet, as 
suggested bv I Shkolnik et al.l (12008 *). Given the short orbital pe- 
riod of the planet, a different approach should be used to search 
for signatures of a possible star-planet int eraction, as in the case 
of CoRoT-2 (cf., e.g., .Pagano et al.ll2009l) . 
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Fig.3. L'pper pane/.- The ME-regularized composite best fit to the out-of-transit light curve of HD 189733 obtained for Q - 0.0. The 
flux is in relative units, i.e., measured with respect to the maximum observed flux along the light curve. Lower panel: The residuals 
from the composite best fit versus time. 



5.3. Activity-induced RV variations 

To simulate the apparent RV changes induced by the distribu- 
tion of active regions derived from our light curve modelling, 
we considered a spectral line with a rest wavelength of 600 nm 
that is close to the isophotal wavelength of MOST observations 
for which our contrast coefficients are given. The spectral reso- 
lution of the line profile is A/AA - 75,000, i.e., comparable to 
that of SOPHIE. 

We fit the observed RV residuals by adjusting an offset value 
and the parameter (Cd) that measures the average reduction of 
the line depth inside starspots. We exclude from the fit the four 
measurements between HJD 2454308.0 and 2454310.0 because 
those residuals show a very steep variation with a decrease of 
the RV of ~ 25 m s"' in only two days, that is impossible to 
reproduce with our model. This fast change and other similar 
rapid variations in the observed RV time series exceed the ac- 
curacy of the SOPHIE velocimetry (cf. Sect. 12.2b . They could 
be due to photospheric velocity fields of several km s ' local- 
ized in active regions and displaying short-term variations, such 
as those observed in the Sun during flares or the emergence 
of ne w magnetic flux, or in syphon flows (e.g., Rtie di et al. I 
1992h. or in the time-v arying features of the Evershed flows (e.g.. 



Cabrera Solana et alTI ^007). We looked at the variabiUty of 



the residuals in the other published RV datasets (iBouchv et al. I 
[2005 ; Winn et al. 2006). They also show RV changes on short 
time scales with night-to-night variations of several tens of 
m s~\ in agreement with the present SOPHIE observations (see. 



e.g., the bottom panel of Fig. 1 in iBouchv et al. 1120051) . The 
chromospheric activity, as monitored by the core emission of 
the Ca II H&K lines, shows remarkable variations on time scales 
as short as one hour, suggesting frequent fl aring events that can 
contribute t o the short-term RV variations dMoutou et alj|2007l : 
Fares et a O20I0h . 

The best fit for the spot model with Q = 0.0 is found for a 
RV offset of 21.6 ± 1. m s"' and (Cd) = 0.25 + 0.05. We plot 
the observed and simulated RV residuals vs. time in the upper 
panel of Fig.|6]together with their differences in the lower panel. 
The eiTor bars of the synthesized residuals take into account the 
effects of the photometric accuracy, the differential rotation, and 
the photometric residuals. The nominal accuracy of MOST pho- 
tometry is ~ 1 . 1 X 10"^ in relative flux units that implies an error 
in the estimated spot area of the order of ~ 2.5 x 10""* lead- 
ing to a RV error of ~ 0.75 m s ' when Eq. (1) in Deso rt et aP 
(|2007|) is applied to compute the maximum deviation. The differ- 
ential rotation changes the longitudes of the starspots during the 
time intervals considered for the photometric modelling. Since 
our modelling assumes a fixed pattern of spots during each in- 
terval of 7.625 days, our synthesized RV variations have sys- 
tematic errors produced by the migration of the active regions 
during those time intervals. We estimate a maximum deviation 
of 1 .5 - 2.0 m s"' depending on the spot migration rate. Finally, 
the small active regions that evolve on timescales of 1 -2 days as 
revealed by the residuals of the photometric best fit (see Fig. |3] 
lower panel) also contribute to the uncertainty of our synthe- 
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Fig. 4. The distributions of the spotted area vs. longitude at the 
labelled times (HJD - 2454000.0) for Q = 0.0. The plots have 
been vertically shifted to show the migration of individual spots 
(relative maxima of the distributions) versus time. The vertical 
short-dashed lines mark longitudes 0° and 360°, beyond which 
the distributions have been repeated to help following the spot 
migration. The long-dashed lines labelled with capital letters 
from A to E trace the migration of the most conspicuous spots 
detected in the plots, respectively (see the text). 



and systematically above for epochs later than 2454310.0. The 
deviations are of about 10 - 15 m s~' for the spot model with 
2 = 5, and are already significant for Q - 1-2. We con- 
clude that the observed RV fluctuations induced by stellar ac- 
tivity favour a spot model with a negligible facular contribution, 
supporting the result found in Sect. |4]for the facular-to-spotted 
area ratio considering only the photometric observations. 

iBoisse et al. I (l201ll) noticed that most of the power of the 
RV variations induced by stellar activity is concentrated at the 
rotation frequency and its first two or three harmonics. Although 
the number of data points is limited to 23 in our case, we applied 
their technique to test whether the subtraction of the synthesized 
RV variations is capable of reducing the impact of stellar activ- 
ity on the RV measurements. We plot in Fig.|2]the Lomb-Scargle 
periodogram of the time series of the RV measurements simulta- 
neous with MOST photometry (solid line) and the periodograms 
of the synthetic RV time series (dash-dotted line, upper panel) 
and of the residuals obtained after subtracting the synthetic se- 
ries from the observations (dotted line, middle panel), exclud- 
ing in all the cases the four points between HJD 2454308.0 and 
2454310.0. All the periodograms have the same normalization 
to allow an immediate comparison of their power levels. The ob- 
served RV variations show the maximum power close to the ro- 
tation frequency with the successive harmonics having a slightly 
increasing power up to the third. This effect is likely due to the 
presence of multiple, partially unresolved peaks among which 
the power is split, as clearly evident in the case of the first har- 
monic. Therefore, the peak of the third harmonic is higher than 
those of the first and the second because it is narrower The spec- 
tral window of the time series is plotted in the lower panel of 
Fig. |7] and shows the presence of significant sidelobes, notably 
at 0.07 d ' from the main peak. This, in combination with the 
remarkable differential rotation of HD 189733, can explain the 
multi-peaked features in the power spectrum, particularly evi- 
dent in the signals of the first and second harmonics. 

When we subtract the synthesized RV variations, the power 
is significantly reduced at the rotation frequency and its first and 
second harmonics, although some residual power is still present 
between the fundamental frequency and the first harmonic and 
between the first and the second harmonics, probably as a conse- 
quence of the multiple frequency nature of those peaks in com- 
bination with the spectral window. The reduction at the third 



sized RV values. Since their estimated areas are of the order of 
0.2 - 0.3 percent of the stellar disc, the estimated fluctuations 
are of ~ 1.8 m s"' using the model of Sect. 13.31 In conclusion, 
the overall standard deviation of our synthesized RV values is 
-2.8 -3.0 ms-'. 

The agreement between our synthesized RV variations for 
Q - and the observations is generally good, although our 
model cannot reproduce the oscillations shown by the data on 
timescales as short as 1 - 2 days because it assumes a fixed 
configuration of active regions that is changed with a timestep 
of 3.8 days. Assuming the standard deviation of the RV mea- 
surements given by the SOPHIE reduction pipeline, we find a 
reduced x^ - ^ -57 for the best fit of the 23 simultaneous data 
points. On the other hand, considering the observed RV fluctu- 
ations on timescales shorter than 1 .0 day, we derive a standard 
deviation of 5.55 m s"' giving a reduced ;^'^ - 1.18. 

Including a facular contribution, i.e., considering models 
with 2 > 0, leads to an increase of the reduced ;i'^ of the best fit. 
This happens because faculae produce an increase of the ampli- 
tude of the RV variations in our model leading to values that are 
systematically below the observations before HJD 2454308.0 
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Fig. 5. The total spotted area as derived from the regularized ME 
models vs. time for Q - 0.0. 
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harmonic after the subtraction of the synthetic RV time series 
is less remarkable. This is expected because our approach can- 
not model the variations on timescales shorter than 4-5 days. 
Unfortunately, the small number of simultaneous RV and photo- 
metric observations limits our capability to test the performance 
of our approach. Nevertheless, the present application demon- 
strates that it can be highly advantageous even for a small sam- 
ple of data with limited precision, i.e., a photometric precision 
of a few percent of the flux modulation and a RV precision of 
3 - 5 m s"'. 

6. Discussion and conclusions 

We have applied th e M E spot model l ing me thod introduced by 
iLanza et all ( l2009al) and lLanza et"an ( 12009 b') to the MOST pho- 
tometry of HD 189733. It is among the brightest and most ac- 
tive stars accompanied by a transiting hot Jupiter and displays 
activity-induced RV variations with an amplitude reaching sev- 
eral tens of m s"' . The simultaneous RV measurements collected 
with SOPHIE offer us a unique opportunity to compare the mod- 
els proposed to account for the RV variations induced by activ- 
ity with the observations. Both the MOST photometry and the 
SOPHIE RV measurements have errors comparable with those 
expected with larger telescopes on fainter stars. In the case of 
CoRoT or Kepler, it is possible to reach a photometric relative 
accuracy of 1.5 x lO"'* and 2.0 x 10"^ on a G2V star of mag- 
nitude Y - 12 in one hour integration time, respectively. Using 
HARPS, we can improve the radial velocity accuracy down to 
2 - 5 m s"' on stars of that apparent magnitude, making our test 
particularly relevant for those targets, especially for the most ac- 
tive ones, given that HD 189733 is among the most active plan- 
etary hosts found so far. 

Our spot models show several active regions with lifetimes 
comparable or longer than the duration of MOST observations, 
i.e., ~ 30 days, that evolve remarkably and migrate at different 
rates with respect to a reference frame rotating with the mean 
stellar rotation period. This is interpreted as evidence of sur- 
face dififerential rotation with a relative amplitude of AQ/Q - 
0.23 ± 0. 10, in agreem ent with t he result o btained from Doppler 
Imaging dMoutou et al. 2007: Fa res et alTl KOlO). It is interesting 
to note that, although the MOST observations cover only ~ 2.5 
stellar rotations, the differential rotation signal is fairly evident 
thanks to its large amplitude. Moreover, a population of small 
active regions with a surface of 0.1 - 0.3 percent of the stellar 
disc and typical lifetimes of 1 - 2 days is observed, as indicated 
by the oscillations in the residuals of the photometric model (cf. 
Fig.©. 

We introduce a method to synthesize the RV variations in- 
duced by stellar activity from the longitude distribution of the 
stellar active regions derived by the ME spot modelling. The 
model assumes that the stellar active regions have the same con- 
trasts and turbulent convection properties of the solar active re- 
gions. Starspots are assumed to be close to the stellar equator 
because we cannot derive their latitudes from the high-precision 
photometry. Indeed, in the case of HD 189733 the spot latitudes 
could be estimated from the migration of the active longitudes, 
given the differential rotation law of Fares et al. (2010), but the 
results are uncertain and we prefer to compute the RV variations 
adopting nearly equatorial spots to test our approach in the gen- 
eral case when the latitudinal dependence of the rotation period 
is unknown. 

Our model has only two free parameters, i.e., an offset ap- 
plied to match the zero points of the synthesized and observed 
RV scales, and a coefficient (Cd) measuring the average reduc- 




Fig. 7. Upper panel: Lomb-Scargleperiodogramof the observed 
RV variations produced by stellar activity excluding the four 
data points between HJD 2454308.0 and 2454310.0 (solid line) 
and the periodogram of the coiTesponding synthetic time series 
(dash-dotted line); Middle panel: Lomb-Scargle periodogram of 
the observed RV variations produced by stellar activity exclud- 
ing the four data points between HJD 2454308.0 and 2454310.0 
(solid line) and the periodogram obstained after subtracting the 
synthesized RV variations from the observed RV time series 
(dotted line). The vertical dashed lines mark the rotation fre- 
quency and its harmonics for a rotation period of 11.953 days. 
Lower panel: The spectral window of the time series of the RV 
observations considered in the upper panel. 



tion of the depth of the spectral lines in dark spots. They are 
derived by a best fit to the RV observations. The synthesized 
RV variations obtained from the starspot distributions show a 
remarkably good agreement with the simultaneously observed 
RV variations, although the latter display remarkable changes 
on timescales as short as 1 - 2 days that cannot be accounted for 
by our model whose spot distributions can be updated only with 
a time interval of ~ 4 days. This is a limitation of every model 
based on the rotational modulation of spot visibility because it 
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Fig. 6. Upper panel: The observed RV variations due to stellar activity (filled dots) and the corresponding synthesized variations 
(open diamonds) versus time for the spot models with 2 = 0. The vertical dotted lines mark the time interval with the steepest 
variations that cannot be fitted by our model (see the text). Lower panel: The difference between observed and synthesized variations 
vs. time. 



requires a minimum time interval, not too short in comparison 
with the stellar rotation period, to map the brightness inhomo- 
geneities on the stellar surface. 

The RV variations set also a constraint on the facular effects 
suggesting that dark spots are dominating both the photomet- 
ric and the activity-induced RV variations, at least during the 
time interval covered by the present observations. Of course, this 
does not exclude the possible presence of bright faculae in the 
photosphere of HD 189733, but their contribution cannot be de- 
tected in the available data. We show that the best value of Q, 
the facular-to-spotted area ratio, is Q = 0, with an acceptable 
range extending from ~ to ~ 5 when only the modelling of the 
optical flux variations is considered. Considering the modelling 
of the RV variations, an upper limit of Q ~ 2 - 3 is derived. In 
the case of the Sun, the best value of 2 is 9 dLanza et al. 2007|). 
Thus our results indicate a lower relative contribution of the fac- 
ulae to the light variation of HD 189733. The amplitude of the 
rotational modulation of our star was ~ 0.015 mag during the 
present MOST observations and ~ 0.03 mag during the 2006 ob- 
servations described in Miller-Ricci et al. (2008), i.e., ~ 5 - 10 
times that of the Sun at the maximum of the eleven-year cycle. 
This indicates that HD 189733 is remarkably more active than 
the Sun, which may account for the undete ctable facular contri- 
bution to its light variations, as suggested bv lRadick et al.l (Il998l) 
and Lockwood et al. (2007). 

We conclude that the approach presented in this paper can 
reduce the power of the activity-induced RV variations at the ro- 



tation frequency and its low-order harmonics by a factor ranging 
from 2 to 10, which represents an important improvement for the 
detection of exoplanets corotating with the star or with an orbital 
frequency in the 2 : 1 ratio with the rotation fr eque ncy, as sug- 
gested by the statistical studies of iLanza I (1201 Ol) and lLanza et al.l 
(201 1). Thanks to the high-precision space-borne photometry al- 
ready available with MOST, CoRoT, and Kepler, it is possible to 
map the longitudinal distribution of active regions in late-type 
stars and apply the method presented in this paper to reduce re- 
markably the impact of stellar activity on their RV jitter allowing 
us to confirm the detection of telluric planets or refine the mea- 
surement of their mass (see, e.g., the application to CoRoT-7 in 
ILanza etal.ll2010h . 
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